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Abstract 

We reconsider the problem of diffusion of particles at constant speed and 
present a generalization of the Telegrapher process to higher dimensional 
stochastic media (d > 1), where the particle can move along 2 d directions. We 
derive the equations for the probability density function using the "formulae 
of differentiation" of Shapiro and Loginov. The model is an advancement over 
similiar models of photon migration in multiply scattering media in that it 
results in a true diffusion at constant speed in the limit of large dimensions. 

I. INTRODUCTION 

Diffusion theory is an old subject and has been applied with enormous success in many 
fields of Physics, Chemistry and Biology. But pure diffusion of particles is actually unphysical 
as it neglects the inertia of the particles. It is a Wiener process for the spatial position of 
the particles and holds strictly in the limit when the trasnport mean free path /* — > and 
the speed c — > oo such that the diffusion coefficient D = d*/3 is a constant. In other words 
the particle scatters at every point and has no directional memory or persistence. The 
persistence becomes important in many problems when processes at short length and time 
scales become important as in polymer physics and radiative transport in random media. 
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The incoherent energy transport of light in random media described by the radiative transfer 
equation [[T],|2| can be approximated as a diffusion process in a coarse grained description. 
But again, the diffusion approximation fails here at short length- and time-scales ||, not 
because of any finite photon mass, but due to the anisotropic (forward) scattering nature 
of the finite size scatterers. This becomes particularly important in the biomedical imaging 
and diagnostic applications with NIR light @] when the ballistic and 'snake' (near ballistic) 
light in highly forward scattering media are used. Since the general analytic solutions to the 
radiative transfer equation that take into account the full scattering properties are not known 
even for simple geometries, it becomes important to develop simple models that incorporate 
the physical persistence of the problem into the stochastic process. Further in the case 
of light, assuming the distance between the scatterers is much larger than the wavelength 
of the light, the speed of the photon is constant in between the scattering events. While 
approximately describing light as a particle undergoing stochastic motion, this constancy of 
the speed should be preserved at the very least. 

Recently we successfully described this diffusion at a constant speed of photons as a 
diffusion on the velocity sphere [Q. Among the other models proposed to deal with the 
intermediate range of length- and time-scales between the ballistic motion and diffusive 
transport, the attempts to generalize the Telegrapher equation are important. The Telegra- 
pher equation is exact in one dimension and describes the diffusion of a particle whose speed 
is fixed, i.e., the velocity can take on only two values ±c 0. The Telegrapher equation, i.e, 

d 2 P dP 2 d 2 P 

where P is the probability distribution function and T/2 is the mean scattering rate, is 
a combination of the wave equation describing the inertial (persistent) aspect and of the 
diffusion equation describing the stochastic aspect. This equation has found wide application 
in many fields and was first considered by J.C. Maxwell |||9[], more than a century ago in 
his attempt to describe heat conduction from basic kinetic theory. He discarded the inertial 
term, however, on the grounds that it would be important only at extremely short time 



scales. It also has since been shown to describe the 'second sound' in liquid Helium II, and 
classical mesoscopic diffusion obeying the Maxwell-Cattaneo law instead of the usual Fick's 



law 1C 



Following a suggestion of Ishimaru 11 1, a generalization of the Telegrapher equation to 



higher dimensions in a heuristic manner was attempted [T2| by simply replacing d 2 P/dx 2 in 
Eqn.(7.1) by V 2 P, to describe photon migration at short length scales, by including some 
ballistic aspects. This ad-hoc generalization appeared to be quite successful to describing 
photon migration as it preserved causality and did much better than the diffusion approx- 



imation to describe photon transport in absorbing media [|i2yi3| . This was also extended 



to studies of Diffusing Wave Spectroscopy in thin samples |14| . It was, however, shown by 
comparing it with Monte-Carlo simulations that this generalization furnished no better an 
approximation than the diffusion approximation in higher dimensions ||15|| . The in-principle 



weakness of this 'ad- hoc' approach was demonstrated by the fact that the photon probability 
density evolving under this equation becomes negative in two dimensions for the simplest 
case of an unbounded (infinite) medium at short times (t ~ t*) when the ballistic aspects 
of transport are most important. In fact, the negativity of the solution to this equation is a 
generic property in even- dimensional spaces |T5|JiT)[ . 



In this paper, we reconsider the problem of diffusion of photons at constant speed and 
present a generalization of the Telegrapher process to higher dimensional turbid media (d > 
1), where the photon can move along 2 d directions along the diagonals of a (^-dimensional 
hypercube. We derive the equation for the probability density function using the "formulae 
of differentiation" of Shapiro and Loginov ||17|| , by considering a correlated random walk 



at constant speed. We show that a partial differential equation of order 2 d results for the 
probability distribution function in ci-dimensions. Our model is an advancement over the 



earlier models of Boguna et al. [18|,|T9|], where the photon could only move along the 2d 
directions along the axes, and results in a true diffusion at constant speed in the limit of 
large dimensions. Our work brings out certain features that were not recognized in earlier 
work. Light in the stochastic medium is considered to be a particle on which the medium 



exerts fluctuating forces. Each scattering event only changes the direction of the photon 
without affecting the speed of propagation. 



II. THE TELEGRAPHER PROCESS IN ONE DIMENSION 

Let us first consider the dynamics of a particle executing a random walk in one dimension 
while moving with constant speed c. This would describe the motion of light in a disordered 
fibre, or of electrons on the Fermi points in a one-dimensional disordered wire, if we neglect 
the wave nature and the consequent Anderson localization (Strictly speaking, this description 
would not hold for 1-D where all the quantum states are localized states. For sample lengths 
much smaller than the localization length, however, the transport is almost diffusive). The 
velocity v(t) of the particle is a random function of time such that it can take only two 
values ±c, i.e., a Dichotomic Markov process. If T/2 be the transition probability per unit 
time between these two values of the velocity (( ) indicate averaging over the disorder), 

<«(*)> = 0, (2) 
(v(t)v(t')) = c 2 exp(-r|t-t'|), (3) 

i.e. the velocity is exponentially correlated in time. We note that the stochastic Langevin 
equation for the displacement x — v(t) gives: 

(x) = 0, (4) 
(* 2 ) = 2 -f(t- 1 -^f^), (5) 

i.e., the behaviour at long times (t — > oo) or very large scattering strengths (large T) is 
diffusive ((x 2 ) ~ t) and at short times (t — > 0), the behaviour is ballistic ((x 2 ) ~ t 2 ). 

Next we will derive the equation for the probability distribution function. Let H(x; t) 
be the phase space density of points in the x-t phase space. Now, II satisfies the Stochastic 
Liouville equation : 

dU d 



Averaging over all the realizations of the random function v(t), by the van Kampen lemma 
|2(|, the probability distribution P(x;t) = (Tl(x;t)). We also define W(x;t) = (v(t)IL(x;t)) 
and obtain, 

dP dW 



dt dx 



0. (7) 



Now using the "formulae of differentiation" of Shapiro and Loginov p7| , for a Dichotomic 
Markov process, 

Eliminating W from the above equations, we obtain 

d 2 P r dP__ 2 &P_ 

i.e. the Telegrapher equation for the probability distribution function P. This is an exact 
description of the motion of a particle with constant speed in 1-D. 

The solutions to Equation. @ are well known and, in an infinite medium given by |]6|JT2 



p( X)t ; x = 0,t = 0) = [—) e- lt ' 2 [I {y) + 2Tt-^] 6(ct - \x\), (10) 

where y = (T /2c)v / c 2 t 2 — x 2 , Iq and I\ are the modified Bessel functions of order zero and 
one respectively, and 6 is the Heaviside step function. Note that the solution is zero for 
| x | > ct and thus causality is preserved. The solution indicates that the particles spread 
out symmetrically from the origin, half of them to the left and the other half to the right, 
with a front velocity of c beyond which there are no more particles and heap up at the 'light 
fronts' where there are N/2 exp(— Tt/2) particles (if there are iV particles altogether). 



III. GENERALIZATION TO HIGHER DIMENSIONS 

Now, we seek a generalization of the Telegrapher process process to higher dimensions. 
The simplest way of doing this is to make every orthogonal component a dichotomic Markov 
process which take the values ±c. Thus, in d dimensions, the particle is seen to move along 
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the diagonals of the d dimensional hypercube with a speed ydc. Here the real space is 
continuous, while the velocity space is discrete and can take only 2 d discrete values in d 
dimensions. This is the model phase space that we consider the photon to execute a random 
walk in. For example, in Fig. [I], we show a possible trajectory in the real space for d = 2 
and a particle released from the origin at t — 0. 
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A. The equations for the generalized Telegrapher process 



In two dimensions, we will consider both the x and y components of the velocity of the 
particle to be independent dichotomic Markov processes, 

v x (t) = cxi(t), (11) 

v y (T) = c X 2(t), (12) 

(Xi(t)) = 0, (13) 

{Xi(t)xs(f)) = Si S e X p(-r\t-lf\), (14) 

where Xi(t) are unimodular processes. The particle is thus seen to move along the four 
directions (±i, ±j) with a constant speed y/2c. Now using the stochastic Louiville equation, 
We define the averages P(x,y;t) = (U(x,y;t)),W x = (v x (t)U(x,y;t)),W y = (v y (t)U(x : y;t)) 
and W xy = (v x (t)v y (t)U(x : y;t)). We generalize the "formula of differentiation" of Shapiro 
and Loginov to n-independent dichotomic Markov processes as 



— ( V!V 2 - • •V n U.[v 1 ,V 2 ,- • •V n \)v 1 ,v 2 ,-,vn 



d 

dt 

= -(r\ + r 2 h h r n )(w!t; 2 • • • u n n[ui, v 2 , ■ ■ ■ ,v n ]) vl:V2: ... Vn 



+ [VlV 2 ■ ■ -V n — )v 1 ,v 2 ,-,v n , (15) 

where Vi(t) are independent dichotomic Markov processes, Tj are the respective transition 
rates and II is a functional of v±, v 2 , ■ ■ •, v n . Using the above, we obtain the following closed 
set of equations : 

dt dx dy 

(19) 



dWy dW X 



dx dy 

Eliminating W x , W y , W xy from the above set of equations, we obtain for the probability 
distribution function P(x,y;t) : 



d (d 



P = 



(20) 



dt \dt 

By performing a 7r/4 rotation of the space axes and rescaling the speed to c, this equation is 
seen to be the same as the one derived by Boguna et al. fl8| by a different approach. Unlike 
the " ad-hoc Generalized Telegrapher equation" of Durian and Rudnick fll2l , this partial 
differential equation is of fourth order involving all space and time derivatives. 

Similiarly, in three dimensions, we consider the x, y and z components to the velocity 
to be independent dichotomic Markov processes. Again, following the above procedure, we 
obtain the closed set of eight equations : 



dP dW x dW v dw- 



elt dx 



+ 



dy 
dW x 

~dT 

dW, 



+ 



dz 







+ TW X = -c 2 — 



w +rw * 



—c 



dW z 
dt 



+ TW V 



9 -^ + 2TW xy 



dt 
dW, 



dt 
dW z 



dt 



+ 2TW yz 
+ 2TW ZX 



—c 



dW, 



xyz 



ot 



+ 3TW : 



xyz 







(21) 


dP dW xy 
dx dy 


dW zx 
dz 


(22) 


dP dW xy 
dy dx 


dW yz 
dz 


(23) 


dP dW zx 
dz dx 


dW yz 
dy 


(24) 


~dW y | dW x ~ 
dx dy 


dW xyz 
dz 


(25) 


~dW y ! dW z 

dz dy 


dW xyz 
dx 


(26) 


~dW z ! dW x 

dx dz 


dW xyz 
dy 


(27) 


~dW yz | dW zx | dW xy 
dx dy dz 


(28) 



It is possible to obtain a cumbersome-looking partial differential equation for P(x,y,z;t) alone, 
similiar to Equation (|20|) by eliminating the other functions from the above coupled set of 
differential equations. However, no extra information results and it will not be presented 
here. 

The partial differential equation for P(r;t) alone, in general, will be of order 2 d (there 
are 2 d independent first order coupled differential equations for d dimensions). The set of 
coupled first order differential equations (f2"T|-|2"8"|) offer a very convenient factorization of the 
2 d dimensional equation satisfied by P(f; t) and are a more convenient starting point for 
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numerical calculations of the solutions. We note that the above equations are linear with 
constant coefficients and can easily be solved by taking Laplace transforms with respect 
to time and space variables. Inverting the solutions obtained back to the real space-time, 
however, is non-trivial and only a few characteristic quantities such as the moments of the 
residence times have been calculated [01 for a similiar model. 



B. Absorbing boundary conditions 

For the above system of equations, absorbing boundary conditions can be easily and 
rigorously applied for this set of equations. For an absorbing boundary at x — 0, with the 
stochastic medium occupying the negative semi-infinite half-space, the appropriate boundary 
conditions corresponding to no incoming flux are 

W x (x = 0, y, z- 1) = cP(x = 0, y, z; t), (29) 

W xy (x = 0, y, z; t) = cW y (x = 0, y, z; t), (30) 

W zx (x = 0, y, z; t) = W z (x = 0, y, z; t), (31) 

W xyz (x = 0, y, z; t) = cW yz (x = 0, y, z; t), (32) 

and free boundary conditions on the other functions P, W y , W z , W yz . These conditions are 
equivalent to the integral boundary condition d/dt Pf~° P(x,y, z;t)dx = cP(x = 0,y,z;t) 
on P(x, y, z; t) alone. 



C. Projected motion along any axis and angular non-symmetry of the model 

In higher dimensions (d > 1), the ad-hoc generalized Telegrapher equation viz. + 
— c 2 V 2 P = is indeed obtained only if higher order velocity correlations are neglected, 
i.e., terms such as W xy = (v x Il), W xyz etc. are set to zero. This is not correct especially at 
short times, when we expect the velocity components to be correlated to a quite some extent. 
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However, it is easily seen that the marginal probability distribution for the projected motion 
along one of the axis p(xi] t) — J P(x 1 , x 2 , ■ ■ ■ , Xd)dxi dx 2 ■ ■ ■ dxd satisfies the Telegrapher 
equation ^ + — c 2 |^§ = 0. The partial differential equation for P(f; t) alone, in general, 
will be of order 2 d (there are 2 d independent first order coupled differential equations for d 
dimensions) corresponding to 2 d directions. 



There are some subtle differences between our model and that of Boguna et al fl8fl . First 
of all, the number of allowed directions for the photon motion is greater in our model (2 d ) 
than theirs (2d). The reason is that, they consider that the motion of the particle to be along 
the axes, while in our case, the motion is along the diagonals of the d dimensional hypercube. 
They do not obtain a Telegrapher equation for the marginal probability distribution for the 
projected motion in general as we do. We always have a Telegrapher process along any 
one axis. This can be best compared in two dimensions by carrying out a ir/A rotation 
of the axes in Eqns. (16-19) and then looking at the projected motion along the axes. We 
obtain the Telegrapher equation + 2T^ — c 2 |^§ = , i.e., only the diffusion coefficient is 
renormalized. This corresponds to the three step Telegrapher process of moving at constant 
speed to the left or the right with a probability 1/4 and being at rest with a probability 
1/2. In higher dimensions (d > 2), the diagonals of the hypercube are not orthogonal and 
the equation obtained for the projected motion along the diagonals is not a Telegrapher 
equation in our case. 

Thus, it is to noted that in these models without angular symmetry, such a description 
of projected motion is non-unique and depends on the direction of the projected motion. It 
should, however, be pointed out that the angular spacing between these discrete directions, 
given by the ratio of the total solid angle to the number of directions, in our model is 

2 d 2 d T{d/2)d^ \2\n l ' 2 (d/2)) ' 1 } 

where V here is the Gamma function. In the limit of large d, the angular spacing decays 
almost exponentially to zero. Thus, in the limit of large dimensions, this process indeed 
describes a genuine diffusion at constant speed. This result is due to the exponential de- 
ll 



pendence on the dimensionality for the number of allowed directions for the photon velocity 
and is not obtained by Boguna et al. 

IV. KUBO- ANDERSON LIKE STOCHASTIC PROCESSES 

In passing, we would like to touch upon the possible generalization to more complex 
stochastic processes, which demonstrates the power of the current approach. It is a simple 
matter now, to write down the equations for probability distribution function for a Kubo- 
Anderson like process, given by the sum of n-independent dichotomic Markov processes (in 
one -dimension), i.e., v(t) = Vi(t) +v 2 (t) + • — \-v n (t), where Vi(t) can take on the values ±c 
and (vi(t)) = 0; (vi(t)vj(t')) = c 2 exp(— T\t — i'|)<%. The structure of the equations forn = 2 
and n = 3 remain the same as Eqn. (16-19) and Eqn. (21-28) respectively, with only the 
derivatives d/dy and d/dz both replaced by d/dx. It is immediately seen that the case of 
n = 2 corresponds to the 3-step Telegrapher processes described above. The generalization 
to higher n is obvious. In general, one obtains n + 1 coupled partial differential equations 
for the sum of n-independent dichotomic Markov processes. 

V. CONCLUSIONS 

In conclusion, we have developed a particular 'generalization' of the Telegrapher process 
to higher dimensional (d > 1) stochastic media which could be potentially useful for studying 
photon migration in turbid media as it rigourously preserves the photon speed to be constant 
between the scattering events. In comparision to the model presented in Ref. ||, where the 
photon's random walk was modelled as a diffusion on the velocity sphere, we have a model 
phase space here, where the photon can move only along the 2 d directions of the diagonals 
of the d dimensional hypercube. It is admittedly an artificial phase space, but one having 
an appreciable directional persistence unlike the zero-persistence diffusion theory. However, 
the model does not have angular symmetry. In 1-D, there are only two directions and hence, 
the Telegrapher equation is exact. On the other hand, in the limit of very large dimensions, 
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the angular spacing between the directions tends to zero almost exponentially, and again 
the process is indeed a genuine diffusion-at-a-constant-speed. It has been shown that the 
equation for the projected motion along any hypercube axis is a 1-D Telegrapher equation, 
though it is non-invariant under an arbitrary rotation. Further, the ad-hoc generalized 
Telegrapher equation in higher dimensions ]T2| , P^| is recovered when higher order correlations 
are neglected. The power of this approach is demonstrated by deriving the equations for a 
sum of n-independent Markov processes. 
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